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Abstract 

We report the consequences of implementing momentum dependent in- 
teractions (MDI) on multifragmentation in heavy-ion reactions over entire 
collision geometry. The evolution of a single cold nucleus using static soft 
equation of state and soft momentum dependent equation of state demon- 
strates that inclusion of momentum dependence increases the emission of 
free nucleons. However, no heavier fragments are emitted artificially. The 
calculations performed within the framework of quantum molecular dy- 
namics approach suggest that MDI strongly influence the system size 
dependence of fragment production. A comparison with ALADiN ex- 
perimental data justifies the use of momentum dependent interactions in 
heavy-ion collisions. 

1 Introduction 



The heavy-ion collisions have always played a fascinating role in exploring vari- 
ous aspects of nuclear dynamics such as fission- fusion [H El [3] , multifragmenta- 
tion m m m [zi m ei uni en ng, collective now on na cm cm as wen 

as particle production [SI [TS1 QI5] etc. The wider energy spectrum available due 
to modern accelerator technologies has become a powerful tool in experimental 
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and theoretical studies in exploring the nature of hot and compressed nuclear 
medium via heavy-ion collisions. 

It is well accepted that the outcome of a reaction depends not only on the 
density but also on the momentum space [TH] HO] . The momentum dependence 
of the equation of state can be extracted from the real part of the optical po- 
tential. This potential is expected to affect those nucleons from target and 
projectile which possess larger relative momenta. The momentum dependence 
of the nucleon-nucleon (n-n) potential is found to affect drastically the collective 
flow observables [UJ [55] and particle production [U [TH1 [Jj5] . It has been shown 
that observables related to the particle production e.g. 7r, k, A yields, nd/n p 
ratios etc. are strongly influenced by the momentum dependence of the n-n 
interaction [TSl [19] . A strong influence of momentum dependent interactions 
was also observed on fragment flow for Ei a b > 400 MeV/nucleon. In higher in- 
cident energy regime, momentum dependent interactions (MDI) cause stronger 
reduction in the number of n-n collisions leading to more pronounced transverse 
flow. The momentum dependent interactions, therefore, increase the mean free 
path of nucleons, and consequently affect the stopping and thermalization of 
nuclear matter [2TJ [23] . 

On the contrary, very few attempts exist in the literature, that shed light on 
the consequences of implementing momentum dependent interactions in frag- 
mentation [TUJ [52] . One of the basic problem with MDI is the strong repulsion 
created in the nuclear environment. As a result, nuclei propagating with MDI 
tend to be destabilized and start decaying via emission of free nucleons and clus- 
ters quite early during the reaction. Interestingly, the role of momentum depen- 
dent interactions depends crucially on the impact parameter of the reaction. It 
is found to enhance the energy of disappearance of flow in central collisions [15] , 
whereas it reduces the energy of disappearance of flow in peripheral collisions 
[TB] ITT] . Similarly, MDI reduce the production of fragments in central collisions 
whereas it enhances the same in peripheral collisions [22] . However, one always 
remained concerned about the stability of nuclei propagating with momentum 
dependent interactions. Even a use of cooling procedure via Pauli potential is 
also reported in the literature [7J. Our present aim, therefore, is to investigate 
the stability of cold nuclei propagating under the influence of momentum de- 
pendent interactions and to see whether one can study fragmentation with MDI 
or not. An attempt to study the system size effects in the presence of momen- 
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turn dependent forces will also be made. We shall also confront our calculations 
with multifragmentation data of ALADiN group [11] which has a rise and fall 
variation with impact parameter. This study is carried within the framework of 
quantum molecular dynamics (QMD) transport model [H0. The QMD model 
and implementation of momentum dependent potential are described in section 
[2j Our results are presented in section [3] and summarized in section [4] 



2 QMD Model and Momentum Dependent In- 
teractions 

The quantum molecular dynamics (QMD) model [HUH] is an n-body transport 
theory that simulates the heavy- ion (HI) reactions between 30 MeV/nucleon 
and 1 GeV/nucleon on event by event basis. It includes quantum features like 
Pauli blocking, stochastic scattering and particle production. Here each nucleon 
is represented by Wigner distribution function of the form: 

/ 4 (r,p,t) = 1 e-^Ml'/VW)^, (1) 
{nil) 6 

where i"i(t) and Pi(t) define the classical orbit, the center of i th Gaussian wave 
packet in phase space which evolves in time. The centers of these Gaussian wave 
packets propagate according to the classical equations of motion [H [6j. The 
interaction part used in the QMD model consists of local Skyrme interaction, a 
finite range Yukawa term and an effective Coulomb interaction among protons 
i.e. 

ytot _j_ yYuk _|_ -y-Coul ^\ 

with local Skyrme interaction consisting of two- and three-body interactions: 

V sk = M(r< - tj) + t 2 6(r t - Tj )S( ri - r fe ). (3) 

Since QMD model is a n-body theory, therefore, Eq. © can be reduced to a 
density dependent potential in the limit of infinite nuclear matter limit as: 

U Sk = a {£)+p(£) 2 . (4) 

Po \PoJ 

The momentum dependence of nuclear mean field is included in Eq. (J2J) via 
momentum dependent interaction as: 

V MDI = Uln 2 \h (Pi - Pjf + ll S ( Ti - rj ) , (5) 
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with parameters t 4 = 1.57 MeV and t 5 = 5.0 x 1(T 4 (MeV/c)- 2 . This param- 
eterization is deduced from the real part of the proton-nucleus optical potential 
which reproduces the experimental data upto 1 GeV/nucleon [U[M|. I n an in- 
finite nuclear matter limit, generalized n-n potential (Eq.@ and ([5])) leads to 
following density and momentum dependent potential (without Coulomb and 
Yukawa terms) : 

U(p, p)=a(£\+0 (j-J + t 4 ln 2 [t 5 (p 2 + 1)] (£) . (6) 

The parameters a, (3 and 7 in Eq.([S]) have to be re-adjusted in the presence of 
momentum dependent interactions so as to reproduce the ground state proper- 
ties of nuclear matter. The parameters corresponding to soft, hard and their 
momentum dependent versions are labelled as S, H, SM and HM, respectively. 
The constants a and [3 give the proper rms radii and binding energies of nuclei 
across the periodic table. Parameter 7 gives us possibility to examine the com- 
pressibility and hence equation of state. The set of parameters corresponding 
to soft (S), hard (H) and their momentum dependent versions SM and HM, 
respectively, can be found in Ref. [5]. 

3 Results and discussion 
3.1 Stability of cold QMD nuclei 

To address the question of stability of computational nucleus in the presence 
of momentum dependent interactions (MDI), we initialize a single cold projec- 
tile using soft (S) equation of state and soft momentum dependent equation 
of state (SM). Earlier theoretical attempts ranging from giant monopole reso- 
nances 25 to nucleosynthesis of heavy elements in mergers of neutron stars [26] 
could be explained if the equation of state (EoS) is relatively soft than when 
it is stiff. Another study concerning the linear momentum transfer occuring in 
central HI collisions also showed that a soft compressibility modulus is needed 
to explain the experimental data [37J HE] ■ These observations motivated us for 
the choice of comparatively softer EoS. We follow the cluster emission pattern 
and rms radii of few computational nuclei. Figure 1 shows the time evolution 
of cold QMD nuclei, namely 58 Ni and 197 Au initialized with S and SM inter- 
actions. The cluster emission is followed for the time span of 200 fm/c. Here, 
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Figure 1: The time evolution of heaviest fragment (A max ), free nucleons, LCPs 
[2 < A < 4], MMFs* [5 < A < 9] and IMFs* [5 < A < A P /3] (A P being the 
mass of projectile nucleus) emitted from a single cold nucleus oi 58 Ni (left panel) 
and 197 Au (right panel). Results obtained with soft (S) equation of state are 
represented by solid lines whereas dashed lines show results with soft momentum 
dependent (SM) equation of state. '*' indicates that heaviest fragment has been 
excluded. 



A max denotes the size of residual nucleus. This should be close to that of par- 
ent nucleus if there is no destabilization of the nucleus. The sizes of parent 
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5S Ni and 197 Au nuclei reduce with the inclusion of momentum dependent in- 
teractions compared to nuclei propagating with static soft interactions alone. 
The SM interactions caused an enhanced emission of free nucleons and light 
charged particles LCPs [2 < A < 4]. However, medium mass fragments MMFs 
[5 < A < 9] and IMFs {[5 < A < A P /S\; A P being the mass of projectile} 
are almost insensitive towards momentum dependent interactions. Superscript 
(*) indicates that heaviest fragment has been excluded from the multiplicities 
of MMFs and IMFs. Only a small fraction is emitted as intermediate mass 
fragments (IMFs). It seems that nucleons close in space are emitted in bulk, 
therefore, leading to an enhanced emission of light clusters. On the contrary, 
very few nucleons, LCPs and heavier clusters are emitted when propagating 
with soft EoS. The enhanced evaporation with MDI is also due to repulsive na- 
ture of these interactions. Does this enhanced emission prohibit one to use MDI 
for fragmentation ? If one sees carefully, majority of mass that leaves the gold 
nucleus (for example, with MDI about 19 units are emitted and (A max ) is close 
to 177) is in the form of free nucleons. In the above gold nucleus, out of 19 units 
about 15 are in terms of free nucleons. In other words, we see that nucleons 
from the surface are emitted and there is no contribution towards the emission 
of intermediate mass fragments. One sees that even with MDI, only 0.15 IMFs 
are emitted on the average. Realizing that as many as 10-12 IMFs can be seen 
emitted in Au+Au reaction [32j . this number with MDI is negligible. 

A survey of the time evolution of rms radii of a single QMD nucleus also 
depicts the same picture. We show in Fig. 2, the time evolution of rms radii 
of 5S Ni and 197 Au nuclei till 200 fm/c. The rms radius of nucleus with SM 
interactions increases gradually compared to that initialized with static soft 
interactions. This behavior reflects that MDI create additional repulsions among 
nucleons which leads to enhanced emission of free nucleons. The rms radii of 
gold and nickel nuclei in soft case shows negligible deviation for the characteristic 
time of HI collision. As discussed above, this enhanced radius is due to the 
emission of free nucleons and not due to the IMFs. Therefore, one can study 
the fragmentation with MDI since the structure of IMFs is not altered by the 
inclusion of MDI. 
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Figure 2: The time variation of rms radii of single cold nuclei of 58 Ni (top panel) 
and 197 Au (bottom panel) using soft (S) equation of state (solid lines) and soft 
momentum dependent (SM) interactions (dashed lines). 



3.2 Heavy-ion collisions and system size effects 

After investigating the behavior of cold nuclei initialized with momentum de- 
pendent interactions, let us study the effect of momentum dependent forces in 
heavy-ion reactions. One of the observables linked with the compression and 
expansion of nuclear matter is the density of fragmenting system. The total 
nuclear matter density is obtained as : 

P(M)=T'p^e-<«.«)-/-. (7) 

Here At and Ap stand for the target and projectile masses, respectively. In our 
approach, average nuclear matter density (p/p ) is calculated in a sphere of 2 
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Figure 3: The average nucleonic density (p/p ) calculated in a central sphere of 
2 fm radius versus reaction time for the central collisions of 58 Ni + 58 Ni (top 
panel) and 197 Au + 197 Au (bottom panel). The results obtained with soft (S) 
and soft momentum dependent (SM) interactions are compared at 50 AMeV 
(left) and 400 AMeV (right). 



fm radius. 

In Fig. 3, we display the time evolution of average nucleon density (p/p ) 
reached in the central region for the head-on collisions of 58 Ni + 58 Ni and 
197 Au + 197 Au at incident energies of 50 and 400 AMeV. The maximal average 
density tends to reduce with inclusion of momentum dependent interactions. 
This happens due to additional n-n repulsions created in the system that pro- 
hibits compression of nuclear matter to a significant level. This difference in the 
behavior of {p/p ) calculated using S and SM interactions diminishes at higher 
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Figure 4: The rapidity distribution dNf rag /dy of free nucleons and LCPs [2 < 
A < 4] as a function of scaled rapidity Y c , m jY pro j/ } Y pro j, being the rapidity of 
projectile for the head-on collisions at 400 AMeV. 



incident energies (400 AMeV). This is due to the fact that in central collisions 
at 400 AMeV, most of the initial n-n correlations are already destroyed and 
matter is already scattered, and therefore, repulsion generated due to MDI does 
not play any significant role. As a result, we do not see much difference in 
average central density reached at higher incident energy. Another important 
quantity related with the initial compression of nuclear matter is the rate of 
binary collisions. We have also checked the collision rate for these reactions 
and its behavior is found to be consistent with the density profile. The rapidity 
distribution of nucleons is another useful tool to characterize the stopping and 
thermalization of the nuclear matter. We have displayed in Fig. 4, the fragment 
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rapidity distribution dNf rag /dy of free nucleons and LCPs for central collisions 
of six different symmetric systems at 400 AMeV. The results are displayed here 
using soft EoS (left panel) and soft EoS including MDI (right panel). The ra- 
pidity distribution is more 'isotropic' and nearly full stopping is achieved in 
heavier systems like Au+Au and Er+Er. In lighter systems, on other hand, 
a larger fraction of particles is concentrated near target and projectile rapidi- 
ties resulting into broad Gaussian shape. This feature can be seen in both S 
and SM cases. The lighter systems, therefore, exhibit larger transparency effect 
i.e. less stopping. Such features are also observed in the experimental data 
of FOPI-group [29]. Based on the experimental observations and theoretical 
trends, one can say that smaller the system, lesser is the stopping. With MDI, 
a slight increase in transparency effect is seen due to lesser stopping of parti- 
cles in longitudinal direction. This happens due to reduction in n-n collisions 
which deflect the fragments in transverse direction. As a results, one obtains 
less particles being stopped in longitudinal direction. The system size effects in 
the production probability of different kinds of fragments has been studied and 
predicted by our group [5]. Here we extend the same study with reference to 
momentum dependent interactions. For this analysis, we simulated the central 
collisions of six symmetric systems 40 Ca + 40 Ca, 58 Ni + 58 Ni, 93 Nb + 93 Nb, 
131 Xe + 131 Xe, 16& Er + 168 Er and 197 Au + 197 Au at incident energies of 50 
and 400 AMeV. We also parameterized the multiplicities as a function of to- 
tal mass of the composite system using a power law of the form: cAJ ot ; A tot 
being the total mass of the system. Figure 5 displays the 'reduced multiplic- 
ity' i.e. multiplicity per nucleon of various kinds of fragments. It is clear that 
the system size effects are more visible in soft equation of state compared to 
soft momentum dependent case. A negative slope obtained for the multiplic- 
ity of free nucleons, fragments of mass A=2, and LCPs at 50 AMeV indicates 
their origin from the surface of interacting nuclei. As we move to momentum 
dependent version, additional break up of n-n correlations leads to enhanced 
emission of free nucleons and light charged particles. As a result, multiplicity of 
MMFs* [5 < A < 9] and IMFs* [5 < A < min {A P /3, 65}] (excluding largest 
fragment A max ) gets reduced at 400 AMeV, indicating the vanishing of system 
size effect with MDI. In higher energy regime, cluster production via emission of 
MMFs* and IM Fs* is strongly suppressed in the presence of MDI. It is worth 
mentioning that earlier studies e.g. see Ref. [30], also reported the momentum 
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Figure 5: The final state scaled multiplicity (calculated at 200 fm/c) of free 
nucleons, fragments with mass A=2, LCPs [2 < A < 4], MMFs* [5 < A < 9] 
and IMFs* [5 < A < rain {Ap/3, 65}] as a function of total mass of the system 
Atot- Results shown here are at incident energies of 50 AMeV (l.h.s) and 400 
AMcV (r.h.s). Open circles depict the calculations with soft (S) interaction 
while solid circles are for soft momentum dependent (SM) interactions. '*' 
means that heaviest fragment has been excluded. 
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Figure 6: The mean IMF multiplicity (Nimf) vs the impact parameter b for 
the reaction of 197 Au + 197 Au at 400 AMeV. The QMD calculations (at 300 
fm/c) using soft EoS (open circles) and soft momentum dependent EoS (solid 
circles) are compared with ALADiN experimental data (filled stars). 



dependent potential to be more repulsive for high momentum nucleons. This 
leads to enhanced emission of free nucleons and LCPs. A similar enhancement 
of the nucleons emission and light cluster production was predicted on inclusion 
of momentum dependent effective interactions in the isoscalar nuclear potential 
and symmetry potential [30| 131 j . Contrary to this, with static soft equation of 
state, the production probability of MMFs and IMFs scale with the system size 
as power law: cA T tot with power factor r close to 3/2. Let us now try to confront 
our calculations with experimental data of ALADiN group [11]. The experi- 
mental data is very fascinating because it has been shown that there is a rise 
and fall of multiplicity of intermediate mass fragments with impact parameter 
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[11] . However universality is observed with mass of the system and with incident 
energies exceeding 400 AMeV. In Fig. 6, we display the multiplicity of interme- 
diate mass fragments as a function of impact parameter using soft (S) and soft 
momentum dependent (SM) equations of state. We see that entire spectrum 
is very well reproduced by the momentum dependent interactions. One should 
also keep in the mind that for central impact parameters, different experimental 
groups like FOPI [S], ALADiN QUE] and Miniball [H] differ significantly in 
the multiplicities of IMFs. Overall, we see a clear need of momentum dependent 
interactions in heavy-ion collisions. 

4 Summary 

Summarizing these findings, we here presented a detailed study on the conse- 
quences of employing momentum dependent potential in multifragment-emission. 
Investigation of a single cold nucleus initialized with soft (S) and soft momentum 
dependent (SM) equations of state reveals that momentum dependent interac- 
tions act as destabilizing factor, which results into enhanced emission of free 
nucleons only. However, no change is seen towards artificial emission of IMFs. 
Further, momentum dependent interactions are observed to weaken the system 
size effects studied at 50 and 400 AMeV. A comparison of model calculations 
with ALADiN data on Au+Au reactions favor strongly the use of momentum 
dependent equation of state in heavy-ion collisions. 

This work was supported by a research grant from Department of Science 
and Technology, Government of India vide grant no. SR/S2/HEP-28/2006. 
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